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I. INTRODUCTION 



The concept of multifractality developed in the last decade has proven to be a powerful 
tool for analyzing systems with complex statistics which otherwise appear to be intractable 
It has found direct application in a wide range of fields including turbulence, chaotic 
attractors, Laplacian growth phenomena etc. [||,[||]. Let us give a simple example of a mul- 
tifractal (MF) phenomenon. On a possibly fractal set X C M d of total size R a field <p(r) is 
given at a microscopic scale t. Then normalized moments of this field may have power law 
scaling behavior for £/R — > 0: 

< <f{r) n > I < <p(r) > n ~ (R/£)- Tn . (1) 

Nontrivial multifractal scaling is found if r n 7^ 0. When the moments are averages over the 
sites of X, ip(r) defines a measure on X and rigorous arguments show that the r n are convex 
from below as functions of n 

Here, we generalize an idea of Cates and Witten Q by deriving the MF spectrum in the 
frames of a field theoretical formalism and make use of renormalization group (RG) methods. 
We relate the MF spectrum to the spectrum of scaling dimensions of a family of composite 
operators of Lagrangian 4 field theory. This gives an example of power of field operators 
whose scaling dimensions show the appropriate convexity for a MF spectrum [|^,p|-|T0|] . 

We thus address a special case of a growth process controlled by a Laplacian field. The 
latter may describe a variety of phenomena depending on the interpretation of the field. For 
diffusion limited aggregation this field is given by the concentration of diffusing particles, in 
solidification processes it is given by the temperature field, in dielectric breakdown it is the 
electric potential, in viscous fingers formation it is the pressure |||. In all these processes 
the resulting structure appears to be of fractal nature and is characterized by appropriate 



fractal dimensions ||11|| . The growth and spatial correlations of the structure are governed 
by spectra of multifractal dimensions 0,0] • In general, the boundary conditions determining 
the field will be given on the surface of the growing aggregate itself. It is this dynamic 
coupling that produces the rich structure of the phenomena and seems to make the general 
dynamical problem intractable. 

Here, we study a simpler case where the fractal structure is given and we look for the 
distribution of a Laplacian field p{r) and its higher moments near the surface of the structure 
0. We will follow the diffusion picture, considering the aggregate as an absorbing fractal, 
"the absorber" . The field p(r) gives the concentration of diffusing particles and vanishes on 
the surface of the absorber. More specifically, we consider the Laplacian field p(r) in the 
vicinity of an absorbing polymer, or near the core of a polymer star. In general, we assume 
the ensemble of absorbers to be characterized by either random walk (RW) or self-avoiding 
walk (SAW) statistics. Multifractal scaling is found for the n-moments (p n (r)) of the field 
with respect to these ensembles. 

This formulation of the problem allows us to use the polymer picture and theory devel- 
oped for polymer networks and stars [0)0 an d extended for copolymer stars PJlQ[]. The 
theory is mapped to a Lagrangian 4 field theory with several couplings jl4Hl6| and higher 



order composite operators [|l3| ,|i^|T0| to describe star vertices. 

Our article is organized in the following way. In the next section we present the path 
integral formulation of the Laplace equation and relate it to a polymer representation. The 
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field theoretical representation and renormalization of this polymer model is discussed in 
section III where we discuss the renormalization group (RG) flow and corresponding ex- 
pressions for the exponents r n . We calculate the multifractal spectrum to third order of 
perturbation theory using two complementary approaches: the zero mass renormalization 
with successive e-expansion (see e.g. ||17|| ) and the massive renormalization group approach 
at fixed dimension [18|. For some special cases we reproduce previous results |7j] that were 



obtained in lower order of perturbation theory. In section [TV] we derive the multifractal 
spectra in terms of series expansions. These we present in both of our RG approaches. The 
resulting series are asymptotic. In section [V] we take this into account and obtain numerical 
values only by careful resummation. In section [VT| we discuss our results and conclude the 
present study. 



II. PATH INTEGRAL SOLUTION OF THE LAPLACE EQUATION AND 

POLYMER ABSORBER MODEL 

In this section we show how to describe the diffusion of particles in the presence of a 
absorbing polymer using a "polymer" formalism that represents both the random walks of 
the diffusing particles and the absorber itself in the same way |HJl9||. Let us formulate this 



problem first in terms of diffusion of particles in time. The probability to find a diffusing 
particle at point r\ at time t if it started at point r at time t — can be described by the 
following normalized path integral: 

G°(r 0l r x ,f) = (5(r«(0)-r )tf(rW(t) -n))^,,,). (2) 

The angular brackets in Eq. @ stand for the following average: 

_ /(---)exp(-Ho(rW,t))d{rW} 
<"™>(' w .*>- Jexp(-H (rW,t))d{rW} ' W 

which is performed with the Hamiltonian: 

The integration in Eq. (0) is performed over all paths r^(r) with < r < t. Note that we 
have absorbed the diffusion constant by a re-definition of time. The unit of the dimensionless 
Hamiltonian 7^ is the product of Boltzmann constant and temperature while that of 
time t is the square microscopic length £ 2 . Spatial boundaries may be included in Eq. (0) 
by restricting the path integral to a subspace. G°(r ,ri,t) obeys the following differential 
equation: 

(A+sr£)<*W>.t)-a < 5 > 

Here, d is the dimension of space. In a given volume V a random walker will visit any site 
after some finite time and we may assume a steady state limit for G°(ro, ri, t) for t = oo. In 
this case G°(r ,r 1 ,t) will become independent of r and its limit will define a field p(r): 
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p(rx) = lim — / dr G°(r Q , r x , t). (6) 

This field p(r) then obeys the Laplace equation: 

Ap(r) = 0. (7) 

We introduce boundary conditions in such a way that the field p(r) equals some constant 
Poo at the extremal volume boundaries (or at |r| = oo ) and vanishes on the absorber. The 
absorber itself we describe by a path r^ 2 \s), < s < S 2 . These boundary conditions are 
implemented by an avoidance interaction w 12 punishing any coincidence of the path of 
the RW and the path of the absorber. The correlation function of a random walk in the 
presence of an absorbing path r^(s) with < s < S 2 may then be written as 

G(r ,r 1 ,S 1 ) = (S(r^(0)-r )5(r^(S 1 )-r 1 ) 

x exp { - ^ jf* d Sl ds 2 5(r^( Sl ) - r^(s 2 )}} Ho{ril);Sl) , (8) 

where we have adopted the notation t = Si. 

We are interested in ensemble averaged moments (p™(r + £)) °f the field i n the vicinity 
of the absorber, i.e. with microscopic £. For the RW ensemble the average is performed with 
respect to the Hamiltonian T-Co(r^ 2 \ 5*2), for the SAW ensemble an additional interaction has 
to be included. The moments we calculate as an ensemble average over all configurations of 
the absorbing polymer choosing the site ro on the middle of the polymer PU|. Formally we 
write these moments for £ — > as 

^ „ m+n 

lim (p>o + 0> = lim ^5- / dr a G* mn 

( r 05 r l! • • • 5 r m+n, Si, . . . , S m+n ). (9) 



0=1 



The normalization Z® mQ takes care of the configurations of the absorber, as explained in the 
next section. The correlation function G* is defined as 



ran 



m+n 



G* mn (r , n,..., r m+n , Si,..., S m+n ) = ( J] 5(r» (0) - r )5(r^ (S a ) - r a ) 

a=l 

{Yfi - j- fl S S I 



Here, the absorbing walk is represented by m = 2 paths r 1 ', r^ 2 - 1 , while the remaining n 
paths represent n random walks, as it is shown in the figure |l|. The interaction matrix u a b 
is in this case given by u a b = {0 if a, b < m or a, b > m; W12 else}. 

The limits in Eq. (P) look rather ill defined at first sight, and indeed they should not be 
taken naively. Also the evaluation of the functional integral (|T0| ) is not defined in this bare 
form. Luckily we have at hand the polymer field theory which has dealt with the problems of 
evaluating these formal expressions PJ . We will show below how the theory is mapped to a 
renormalizable O(M) symmetric field theory in terms of which the limits and a perturbative 
expansion of (|T0| ) make sense. For instance the limit |£| — > may be interpreted as a short 
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distance limit denning a composite operator, while the limit S a > m — ► oo, with 5&< m stay 
finite corresponds to a short chain limit derived in [^] . In the frames of the polymer picture 
we may interpret G* mn {r 0l r 1; . . . , r m+n , Si, ... , S m+n ) as the correlation function of m + n 
interacting walks all starting at point r with end points at r 1; . . . , r m+n . These describe 
what is called a polymer star. The normalized partition function of such a star of m + n 



polymer chains with chain lengths parametrized by S a may be written as ||12yi3 |: 

^ „ m+n 

£*mn{S a } = — -= / | J dr a G mn (ro, ri, . . . , r m _|_ n , Si, ... , S m +n). 

J*mri J , 
a=l 



(11) 



The normalization A/" mn is chosen such that -Z* m n{<S'a}|)i ai ,=o = 1 for vanishing interactions 
and the point r$ is arbitrary. We have studied this problem of polymer stars with general 
interaction matrix u a i, in ||9|,[T0l| . Here we will choose 

Mn if a i b < m 

u a b = { u 2 2 if a, b > m (12) 
U12 else 

This allows us to think of the absorbing paths r' 1 ', . . . , r m as being either of a RW (un = 0) 
or of a SAW (uu 7^ 0) ensemble. With m > 2 we also include in our study the moments of 
the diffusion field near to the core of a polymer star. We included U22 only to ease notation, 
in the present context u 22 = 0. 



III. FIELD THEORY AND RENORMALIZATION 

As is well known, the polymer model may be mapped to the limit M = of O(M)- 
symmetrical Lagrangian field theory |23| . We adopt the formalism developed for multicom- 
ponent polymer solutions which allows us to describe both polymers and interacting random 



walks fL6| . Its field theory is given by the following Lagrangian: 

m+n 



-t I I ti I I lr li 

C{<j> a ,»a} = - / d<V + i^Mr)) 2 ) 

m+n „ 

+ 4! E ^ J ^l{r)<Pl{r). (13) 



The <p 2 terms should be read as scalar products of fields <ft a with M component 



% = E(^) 2 - ( 14 ) 



0=1 



The parameter fi a is a chemical potential conjugated to the chain length variables S a in Eq. 
(S). Correlation functions in this theory are defined by averaging with the Lagrangian L: 

((•••)>£= / V[Mr)](---)ex V [-C{<P a ,fi a }}\ M=0 . (15) 
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Here, functional integration j D[(f) a (r)] is defined in such a way that normalization is already 
included: (1)|£ = 1 if all u a>b = 0. The limit M = in Eq. (|15|) can be understood as a 
selection rule for the diagrams which contribute to the perturbation theory expansion and 
can be easily checked in this context to correspond to interacting polymers in the following 
way: The partition function Z* mn defined in Eq. (|j) is mapped to the field theoretical cor- 
relation function Z* mn via Laplace transforms in the chain length variables S a to conjugate 
chemical potentials ("mass variables") /i a : 

poo 

Z*mn{Ha}= / f[ dS^'^ Z* mn {S a } . (16) 

Jo b 

In terms of the above defined Lagrangian field theory Z* mn is given by 

/m+n 
JJ dr a <f) a (r )<f) a (r a ))\ c . (17) 



a=l 



Our interest is in the scaling properties of these functions. Note that by ( |17D these are 
governed by the spectrum of scaling dimensions of the composite operators n^=T 4>a- To 
extract these dimensions we use RG methods |24|,|25|]. Here, we apply the results of our 
previous approaches to the problem of co-polymer stars P,pH|: massless renormalization 



group scheme with successive e-expansion (see e.g. [pj]]) and massive renormalization group 



approach at fixed dimension |Tq| compiled in a pseudo-e expansion [26]. On the basis of cor- 



relation functions it is standard to define vertex functions Eg corresponding to the couplings 
u ab as well as vertex functions T^™" with insertion of composite operators [] a Explicit 
expressions may be found in [JToJ] . We define renormalization and introduce renormalized 
couplings g ab by: 

(18) 

The renormalizing Z-f actors are power series in the renormalized couplings g a b subject to 
the following renormalization conditions: 

d 

Z (j>a (g aa )—T ( ^(u aa (g aa )) = l } (19) 

Z ab {{g ab })Y { ^ bb {u ab ({g ab })) = fi £ g ab - (20) 

The scale parameter [i is equal to the mass at which the massive scheme is evaluated or it 
gives the scale of the external momenta in the massless scheme. 

In order to renormalize the star vertex functions we introduce renormalization factors 

7*mn u... 

k 

(II < 2 )^Z r nZ(^({^})) = ^ ( 21 ) 

a=l 

where 5n<j> a ls the engineering dimension of the composite operator 

5n* a = (m + n)(|-l)+ 4- e. (22) 
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The renormalized couplings g a b and renormalizing Z-factors depend on the scale parameter 
u. This is expressed by the following RG flow equations: 

V-^9ah = Pab({9ab}), (23) 

/^ln^X(fe>}) = VnMU}). (24) 

Our original problem is described by two sets of walks of different species. In this case only 
the three different couplings tin, U\i and U22 in Eq. ( |i~2"|) remain. We will refer to their 
renormalized counterparts as gu, (722, 9\2 = 9<2\- The corresponding functions fin, P22, P12 
define the flow in the space of couplings. This RG flow was discussed in |I5| , |16| . Its fixed 
points are determined by the set of equations: 

AiG?n)=0, 

P22(g* 22 ) = o, 

M9u,9n,9n) = 0- (25) 
In the space of the three couplings one finds 8 fixed points corresponding to the 



absence or presence of inter- and intra- species interaction. The equations for the fixed 
points of the /5-functions were found to have the following nontrivial solutions: /3 a a(9s) = 
and for a ^ b: P 12 (0,0,g* G ) = 0, M^Ati) = °> M0,98,9v) = °> Mfo&gs) = °> 
corresponding to all combinations of interacting and non-interacting chains. 

The phenomenon we address in this article corresponds to the case of a non-vanishing 
interaction between the two species of walks, while one set has no self-interaction. Thus 
we consider the two fixed points which we call G (gu = #22 = 0, g±2 = 9q) and U (gn = 
9*i 922 = 0,(712 = <7(y). The first (G) corresponds to a set of random walks interacting with 
another set of random walks of a second species and thus describes absorption on random 
walk absorbers, the second ( V) corresponds to a set of random walks interacting with a set 
of self-avoiding walks and thus describes absorption on SAW (polymer) absorbers. 

Having m walks of the first species and n walks of second species we define the following 
exponents in the fixed points G, U: 

V°n = Vu^ a (gn = 922 = O,0i 2 = 9a), (26) 
Vmn = Vu<f> a (9u = 9*, 922 = 0, 9l2 = 9*u) , (27) 

which govern the scaling properties of the partition sum (|TT|). 

The scaling may be formulated in terms of the size R of the absorbing walks while the 
RWs of the diffusing particles are taken infinitely long. This corresponds to a short chain 



expansion ||22|| . We have to normalize the partition function by the number of configurations 
of the absorber given by -2* m o and by the nth power of the first moment (see Eq. ([!])). For 
large R on the microscopic scale I the moments of p(ro) at point r$ in the vicinity of the 
core of the star scale like 

< p(ro) n > ~ l7 ( 7 l7 y n (Ryr™ , , 

< p(r ) > n ~ * mn ' * m0 \ * ml ' * m0 J \J ) ' 
where R = S^ <m and the exponents r mn are given as 
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(29) 



Here, 1/ is the correlation length critical exponent of the walks: v = 1/2 for random walks 
and z/ ~ 0.588 for self-avoiding walks at d = 3. For the fixed point (7 we have 7^ = and 
v — 1/2 for all walks. The scaling near a RW star is then described by inserting the value of 
r\ while the scaling near a SAW star is obtained by inserting rj u . In previous work P1,|T0 
we obtained the exponents 7^ n , 7/^ n in third order of perturbation theory. 



IV. MULTIFRACTAL SPECTRUM 

A widely used characterization for the MF spectrum is the so called spectral function 
f m {oi) @. To obtain this function for the absorption process on the center of a star with 
m legs we analytically continue the set of exponents T mn in the variable n and calculate the 
following Legendre transform 

fm{ot m n) = -T mn + na mn + D rn with a mn = + D m . (30) 



Following the standard definition we have included into Eq. (|50D the fractal dimension D rn 
of the absorber. In particular, this gives the maximal value of the spectral function f m (a mn ) 
to be equal to the dimension D m . However the last is ill defined in the case of absorption 
near the core of a polymer star. That is why in the subsequent analysis we will neglect 
D m in the expressions (|30|). Note as well that we have defined the exponents T mn based on 
ensemble averages for the moments < p(r) m >. This definition deviates from the standard 
approach where an average over the support of p(r) is assumed. One has to expect some 
deviations from the behavior of the standard f{o) which we will discuss below [[/J. To 
obtain the expressions for the spectral function we use the perturbation expansions for the 
r) exponents given to third order both in massless and massive renormalization |9],[10| . These 



exponents are available both in terms of e-expansion and pseudo-e-expansion series. The 
first corresponds to collecting perturbation theory terms of the same power of e = 4 — d. 
In the pseudo-e expansion series |26[ each power of the pseudo-e parameter r collects the 



contributions from the dimension- dependent loop integrals of the same order. In the final 
results the limit r = 1 is taken. Starting from the relations for r mn (|29|) and the spectral 
function (p0|) some algebra results in the corresponding expansions for the MF spectra for 
absorption on stars of random walks (RW) and self avoiding walks (SAW): 



a RW (e 



-m(2n - l)e 2 / 8 + rnfimn + Gn((3) - 12n + 3n 2 (31) 
+5-2m-3C(3))£ 3 /i6, 

f* w (e) = -rrmV/8 + mn 2 ( - 6 + 2n + 2m + 3C(3))^ 3 /16, (32) 

aii w (e) = -9m(2n - 1^ 2 / 12 8 + 3m(l68mn + 54n 2 + 157 (33) 

+180nC(3) - 350n - 84m - 90C(3)^ 3 /2048, 

fl AW (e) = -9mnV/128 + 3mn 2 ( - 175 + 36n + 84m + 90C(3)^ 3 /2048, (34) 
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a RW (r 



-me (l + Ann - 2i x - 2nj r 2 /4 + c^ p r 3 , (35) 

/r (r) = W(- 1 + 2i 1 )r 2 /4 + (36) 

«£f» = -9me(l - 2n + Am, - 2^)^/64 + og£?T* (37) 

/^(r) = 1 + 2z x )r 2 /64 + (38) 

Here £(3) ~ 1.202 is the Riemann zeta function, ii, i 2 are the two-loop integrals depending 
on the space dimension d: at d = 3 i\ = 2/3, i 2 = —2/27. The explicit form of the three-loop 
contributions in Eqs. fl35|)-(]38[) is given in the appendix [A]. 

V. RESUMMATION AND RESULTS 

As is well known, the series of type (|3T| ) - (|38"D, as they occur in field theory appear to 
be of asymptotic nature with zero radius of convergence. However, knowing the asymptotic 
behavior of the series as derived from the RG theory we may evaluate these asymptotic series 
by resummation (see e.g. ||27|| ). To this end several procedures are available differing in the 
amount of information that is used to control the convergence. We extract this additional 



information for the case of our Lagrangian (O) from [16, 28 1 . We expect the following 



behavior of the kih order perturbation theory term Ak for any of the above quantities: 

A k ~ k\ k b {-a) k . (39) 

The constant a for the e-expansion of Lagrangian 4 field theory with one coupling was 
derived in [2S|: a = 3/8. For the unsymmetric fixed point U, where two different couplings 



are present the value a = 27/64 has been proposed ||16|| . We assume here that the same 
properties also hold for the pseudo-e expansion in terms of r. With the above information 
at hand one can make use of the Borel summation technique improved by the conformal 



mapping procedure which has served a powerful tool in field theory calculations (see []27| for 
example). We present some details of the resummation procedure in appendix [B]. 

As stated before, we have chosen a definition for the generalized spectral functions f m (a) 
neglecting the fractal dimension D m of the absorber. In this way we ease the comparison of 
the spectral functions / m (a) for different values of m, the number of legs of the absorbing 
star. We have shifted all maxima of f m (a) to the point known for m = 2 with D m=2 = 2 in 
the RW case and D m=2 = 1.71 in the SAW case. 

The results of the resummation of the series for the spectral functions are presented 
in figures |2|-|| Each point marked by a symbol corresponds to the resummation of both 
fm( a mn) an d Oi mn for a given pair (m, n) where we used a half-integer spacing for the values 
of n. Note that the right wings of the curves correspond to negative n < 0. In this region 
reliable resummations were feasible only for sufficiently large m. We have only included 
resummations that were successful in minimizing the deviation between the second and 
third order resummed values as described in the appendix |B[ 

In Fig.|2| we study the effect of different RG and resummation procedures for /^ w (a) in 
d = 3 dimensions. The most simple approach is to directly insert e — 1 into the e-expansion 
and for the r-expansion to use r = 1 and the d = 3-dimensional values for the integrals . As 
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can be seen from the curves, this will work only for the series truncated at second order and 
for small n, i.e. near the maximum of /2(a) at n = 0. In addition we have performed an 
analytical continuation of our series using [2/1] Pade approximants for the series truncated 
at third order. The symmetry of the Pade approximant holds only in the region shown and 
may be an artifact of the method. On the left wing, where it coincides with the resummed 
results the Pade approximant gives a continuation which is compatible with the estimation 
for the minimal a value a m j n = d — 2 ]7|]. The Pade result is a m \ n {e) = 1.333, a m in(T) = 1.017 
for the RW absorber and a min (e) = 1.250, a min (r) = 1.013 for the SAW absorber, which 
is calculated here only from third order perturbation theory. The Pade approximant, while 
already significantly improving the convergence of the results, introduces some apparently 
artificial singularities. Moreover it does not make full use of the known asymptotics for 
the ^-expansions. We have therefore chosen a more sophisticated method of resummation 
that has proven to reproduce reliable data in many field theoretic applications [p5 , p7| . The 



results of these resummations are again shown for both RG approaches. Note that though 
the results obtained for a n and /^(ctn) for a specific value of n differ in both approaches, 
the same curve /2(a) is described with better coincidence for the left wing of the curves, 
corresponding to positive n. 

Figs. present the resummed MF spectra f m {a) of Brownian motion near general 
m-leg polymer stars in d = 3 dimensions. The family of curves f m (a) appear to approach 
some limiting envelope for increasing m in all cases. This behavior is more pronounced in 
the case of Brownian motion near an absorbing SAW star. This provides evidence that 
the MF spectrum catches rather general properties of the phenomena under consideration. 
Here, for the absorption of diffusing particles on a polymer star the spectrum only slightly 
varies with the number of legs m of the star even in the vicinity of the core of the star. Only 
the absorption on an endpoint (m = 1) proves to be an exception. 

The behavior of the maximum of the spectra may also be studied in terms of the series 
expansion. The original position of the maximum is given by its a-coordinate in e-expansion 
in the following form: 

a m fl = r) m ,l ~ V'mfli (40) 

C = m e/S + . . . , (41) 
o%£ = mil - m)e/S + .... (42) 

For m > 1 the position of the SAW maximum is shifted in opposite direction to that of the 
RW maximum. In e-expansion we find for the curvature at the maximum: 

V/m(a) = -<o> (43) 
1//™ W («) = -me 2 /4(1 - e/2(2m - 6 + 3((3)) + . . . , (44) 
l//f AW («) = -9m £ 2 /64 + . . . . (45) 

Here, we use the notations f' m (a) = d/daf m (a m , n )\ n= o and rj' mn = d / dn7] m ^ n \ n=0 and corre- 
spondingly for higher derivatives. As can be seen also in the plots, the radius TZ m ~ 1/ f^(a) 
of the curvature increases with m for both RW and SAW star. Some asymmetry is also 
present in the plots. It may be more explicitly extracted from the series by considering 

/m(«)//m(«) = W'oYKm , (46) 



10 



/r W («)//r W («) = me/12(l - ^(2m - 6 + 3C(3)) + . . . , (47) 
fm AW ( a )/fm AW ( a ) = me/16(l - e{7m/2 - 175/24 + 15/64C(3)) + . . . . (48) 

This shows that the asymmetry at the maximum decreases slightly with m. The plots seem 
to indicate that it approaches some limiting value. 

From the plots we present here, in general one may deduce that the series for the MF 
spectra for diffusion near an absorbing polymer star possess stable resummations and that 
the shape of the resulting curves is robust against the change of the number of legs m of the 
polymer star while a limiting curve seems to be approached with increasing m. 



VI. CONCLUSIONS 

The present work represents extended results on the MF behavior of Brownian motion in 
the vicinity of an absorbing polymer structure. We extend the ideas by Cates and Witten |7| 
to map this problem to a problem of interacting walks. The former authors used a Fixman 
expansion technique to extract the exponents governing the MF scaling. This approach is 
equivalent to a direct renormalization method and unique to dimensional renormalization 
with e-expansion. The Fixman expansion assumes without proof the renormalizability of 
the quantities corresponding to higher moments of the Laplacian field of Brownian motion. 
Here, we map the problem to an 0(M)-symmetric 4 - field theory relating the above quanti- 



ties to corresponding composite operators for which renormalizability has been proven Ul3 
Furthermore, the scaling exponents of these operators have been calculated in our previous 
work ||,[Tn|. There, it has been shown that the resulting spectra of scaling dimensions of 



these operators display the convexity properties which are necessarily found for multifractal 
scaling but unusual for power of field operators in field theory HQ. 

The extensive RG study in PJlCl with three loop results allows us here to consider the 



general case of Brownian motion near the core of a star polymer with m legs. The higher 
order calculation was improved by resummation to give reliable estimates for the families 
of MF spectra describing the multiscaling of Brownian motion near absorbing RW or SAW 
stars. 

Our results have proven to be equally stable under change of the general RG treatment. 
We applied two complementary RG approaches, the dimensional renormalization with suc- 
cessive e-expansion as well as massive renormalization at fixed space dimensionalities. The 
resummation in particular enabled us to extend the region over which the curves for the MF 
spectra coincide in both RG approaches reflecting the stability of the scheme of calculations. 
Our plotted results (figs. 0-§|) indicate some independence of the spectral function on the 
number of legs of the absorbing polymer star. For higher number of legs the spectrum seems 
to approach a limiting curve. The MF spectra calculated here show most of the common 
features shared by spectral functions that describe a variety of MF phenomena. Let us 
note however that unlike the common definition of the underlying scaling exponents based 
on site averages we here rely on ensemble averages for the moments of the Laplacian field 
of Brownian motion. We average over the configurational ensembles of absorbing polymer 
stars. Only for the case of m = 2 legs of the absorbing star, also a site average definition 
could be used. As has been noted also in [0,^] the ensemble average leads to the possibility 
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of negative values of the spectral function (see figs. 0-0) • Furthermore the fractal dimension 
is not defined for the core of a polymer star. 

Experimentally such absorbing systems are realized in diffusion controlled reactions with 
traps or reaction sites attached to polymer chains. This is described by the irreversible 
reaction A + B — > with freely diffusing particles A and traps B attached along the polymer 



chain |fL9|| . The higher nth moments of the field at some special trap might then describe 
the reaction rate for A n + C —>■ 0. If C is located at the core of a polymer star this system 
realizes all aspects of our study. While extensive MC studies exist for many problems in the 
field of DLA we hope that our detailed calculations might initiate also an MC approach to 
the present problem. The current study also allows an extrapolation to d — 2 dimensions. 
Though one should not expect to find results of a pure two-dimensional approach due to 
topological restrictions in d = 2 that are not present in the perturbative e-expansion [|30 . 



While standard in field theoretical studies of critical phenomena, the resummation tech- 
nique as to our knowledge was not applied in the theory of multifractals. We hope that our 
attempt will attract attention for this possibility in the context of other problems that arise 
in the theory of multifractal measures. 



ACKNOWLEDGMENTS 

This work was supported in part by SFB 237 of Deutsche Forschungsgemeinschaft and 
by the Minerva foundation. Yu.H. acknowledges hospitality of the Institut fur Theoretische 
Physik II, Heinrich-Heine-Universitat Diisseldorf, where this work was completed. 

It is our pleasure to acknowledge discussions with Lothar Schafer, Bertrand Duplantier 
and Alexander Olemskoi. We are indebted to Gleb Oshanin for attracting our attention to 



the references [19 . 



a. 



APPENDIX A: THREE-LOOP CONTRIBUTIONS TO THE MULTIFRACTAL 

SPECTRUM 

Here, we collect the three-loop contributions of the expressions for the Holder exponents 
«mn (i5), (0) and the spectral functions f m (a mn ) fl36|) , (|38|) obtained in the pseudo-£- 
expansion scheme. The expressions read: 

DTI/ TJ1/S / Q 9 / A \ 

a 3ioo P = — ^~ \3i 6 m — 3n + 3m — 6mn + 20ii — 6i 4 + 3z 5 — 3i 6 + 3i 7 — 16zi — 46m x (Al) 
—5 + 624m + 3i$m + Yin — Qi^n 2 — 12i 4 mn + 9n 2 ii + 32nii 2 — 6ijn — Qi^mn 

—Qi§mn + Qi§n + l%i^n — Ylmi\ + 2Anmi\ — Qi^rt^j , 

2 

f^ioop = g — \6 — 2n — <6i^n + 9z 4 — 3i 5 m — Qi A m — 3i e m — 3i 7 — 3i 5 + \2m%x (A2) 

+ 16ii 2 — 23i\ + 6nii — 3m + 3i 6 j , 

a*™ = ^ - 86 - 718mi - 198mz! + 36i 6 m + 126i A m + 36i 5 m - 27n 2 + 54m (A3) 

-108mn + 332h + Ai 2 - 108i 4 + 36z 5 - 36i 6 + 45i 7 - 72i 5 n + 72i 6 n - 90i 7 n 
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+496mi 2 + 81n 2 ii — 72i^mn — 81i 4 ri 2 — 252z 4 mn — 12i§mn — 8^1 + \&n%2%\ 
-8m" 2 + 396nm?i - 248i x 2 + 190n + 270i 4 n 

fg£g = __ (95 - i8n - 54i 4 n + 135i 4 - 36i 5 m - 126i 4 m - 362 6 m - 45i 7 - 36i 5 (A4) 

+ 198m?i + 248ii 2 - 359«i + 54nii - 54m + 36i 6 - 4z 2 + 8i 2 *i) • 

The numerical values of the two-loop (ii, i 2 ) and three-loop 13—1*8 integrals at d = 3 equal 



3J: i! = 2/3, i 2 = -2/27, i 3 = -0.0376820725, z 4 = 0.3835760966, z 5 = 0.5194312413, 



z 6 = 1/2, z 7 = 0.1739006107, i$ = -0.0946514319. 

APPENDIX B: RESUMMATION PROCEDURE 

Here, we briefly describe the resummation method for the asymptotic series that we 
applied in our calculations p7^8|, 32|. The starting point is the expansion for the function 



of interest: 

k 

The coefficients A k are supposed to posses the following behavior 

A k = ck b0 {-a) k k\[l + O{l/k)}, k^oo (B2) 

with known values of constants c, 60, a - The property (|B2|) indicates the Borel sumability 
of the series flBj]) . The Borel resummation procedure takes into account the asymptotic 
behavior of the coefficients and maps the asymptotic series to a convergent one with the 
same asymptotic limit. The procedure is as follows. For the series ( |B 1|) we define a Borel- 
Leroy transform /(e) by 

. _ f(j) F j 

/*W = Ef(jT»Ti)" (B3) 

with the Euler Gamma function T(x) and a fit parameter b. Then the initial series may be 
regained from 



POO 

/ dtt b e- t f B {et). (B4) 
Jo 



f ves (e) 

Assuming the behavior of the high order terms ( |B2| ) one concludes that the singularity of the 
transformed series f B {e) closest to the origin is located at the point (— 1/a). Conformally 
mapping the e plane onto a disk of radius 1 while leaving the origin invariant, 

{l + aef^-l 4 w 

w 



(l + ae)V2 + l' a(l-w) 2 ' 

and substituting this into f B (e), and re-expanding in w, we receive a series defined on 
the disk with radius 1 in the w plane. This series is then re-substituted into Eq. (|B4|) . 
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In order to weaken a possible singularity in the w-plane the corresponding expression is 
multiplied by (1— w) a and thus one more parameter a is introduced [32 



In the resummation 
of the e- 



16.28 



procedure the value of a is taken from the known large-order behavior 
expansion series while a was chosen in our calculations as a fit parameter defined by the 
condition of minimal difference between resummed second order and third order results. 
The resummation procedure was seen to be quite insensitive to the parameter b introduced 
by the Borel-Leroy transformation (|B3|). The above procedure was applied to both the e- 
and pseudo-£-expansion series. 
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APPENDIX: FIGURES 




FIG. 1. n random walks that end at point vq on the absorbing polymer. 




FIG. 2. MF spectrum for diffusion near an absorbing RW. Comparison of different approxima- 
tion and resummation schemes. 
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FIG. 3. MF spectra for diffusion near an absorbing RW star (e expansion). 
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FIG. 4. MF spectra for diffusion near an absorbing RW star (r-expansion). 
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FIG. 5. MF spectra for diffusion near an absorbing SAW star (e expansion). 
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FIG. 6. MF spectra for diffusion near an absorbing SAW star (r-expansion). 
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